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Abstract 

Stationary ID Schrodinger equations with polynomial potentials are reduced to explicit 
countable closed systems of exact quantization conditions, which are selfconsistent constraints 
;JJ ■ upon the zeros of zeta-regularized spectral determinants, complementing the usual asymptotic 
(Bohr-Sommerfeld) constraints. (This reduction is currently completed under a certain vanish- 
ing condition.) In particular, the symmetric quartic oscillators are admissible systems, and the 
formalism is tested upon them. Enforcing the exact and asymptotic constraints by suitable itera- 
tive schemes, we numerically observe geometric convergence to the correct eigenvalues/functions 
in some test cases, suggesting that the output of the reduction should define a contractive fixed- 
point problem (at least in some vicinity of the pure q case). 

1 The Schrodinger spectral problem 

We study the 1- dimensional stationary Schrodinger equation with a real polynomial po- 
tential V{q) on the real axis ([Q-IQ provide some directly related references). 
In reduced form, 

- i^"{q) + {V{q) + A)^(g) = 0, V{q) = +g^ + Viq""-' + V2q''-^ + ■■■ + VN-iq. (1) 



We write the collection of coefficients as v = (f i, ■ ■ ■ , vn-i), absorbing any constant term 
into the spectral parameter A. 

The case N = 2 (singular) is excluded [0] . 

1.1 Polarized boundary conditions 

We ffist cast eq.(|l|) into an eigenvalue problem with asymmetric boundary conditions. 
I.e., we restrict the problem to a half-line on which V is confining: [0, +C)o), keeping the 
square-integrability at +oo, then put a Neumann, resp. Dirichlet, boundary condition 
at the finite endpoint q = 0. Both conditions then define self-adjoint operators H~^, 
resp. H~, which have purely discrete spectra respectively denoted {E2k} and {£'2fc+i} 
for k = 0, 1,2, . . .; these admit an asymptotic expansion which is the Bohr-Sommerfeld 
formula reexpanded in energy, 

Y, h-J/NEr^^^ r^27r{k + 1/2), A; ^ cx) in N, with (2) 

i=o,i,2,... 

^= 2 + iV '■ ^^^ growth order; ^/= |,^^,^/dg = ^ r(-)/r(- + -); (3) 

the higher b^^j/^- are polynomials in the Vji {f < j), e.g., &^_i/Ar = —Avi/N. Crucially, 
/i < 1. 

Up to sign, the eigenvalues read as the zeros of the Fredholm determinants 

A^(A)^=^ n (1 + A/i?.); (4) 

7 even 
^ odd 

these infinite products converge to entire functions of order /i in the variable A, also entire 
in the parameters v 0. 

At fixed (A,'u), let ''p\{q) denote a recessive solution of eq.(|l|) (i.e., a solution exponen- 
tially decreasing as g — i> +oo, which is unique up to normalization). Then V'a(O) vanishes 
simultaneously with A^(A), (and likewise for V'^(O) and A^(A)). However, the relations 
connecting them retain messy factors unless each of A and ip is suitably renormalized. 
The procedure for the determinants is the well-known zeta-regularization, but we will also 
introduce a parallel operation upon the wave functions. 

1.2 The spectral (or functional) determinants 

By Mellin-transforming eq.(H), the spectral zeta functions, [|] 

Z±(s)= Y. V (and Z{s) = Z+{s) + Z-{s)), (Res>/i) (5) 

, even 
'^ odd 

are shown to extend meromorphically to all of C, and to be regular at integer points. 



In particular, the point s = (where Z{0) = {2n) ^bo) is used to define determinants 
as zeta-regularized products, 






detH^ = exp[-Z^{0)] and hkewise, D^{X) = det{H^ + X). (6) 

This normahzation, unhke that in eq.(^, commutes with global spectral translations. 
More explicit formulae are 

D±(A)=exp[-Z±'(0)]A±(A); (7) 

log D^{X) = hm { J2 ^og{Ek + X) + l log{EK + X) (8) 

They will also be used for complex spectra, as justified in 0. 

1.3 Natural WKB normalization for recessive solutions 

A recessive solution of eq. (|I|) admits an exact WKB representation, 

^A,„(g) = M(g)-'/'e"^<'o"^'''^^^' where, for q -. +oo, u{q) ~ U{q) = {V{q) + X)'^^. 

(9) 
This normalization awkwardly depends on the base point qo, so we intend to define it for 

go = +oo: 

Mq) = uiqr'/'eC^-^^'^"'' "^^ er°°«(^)'^''^,,o(g). (10) 

However, only u — 11 is integrable at infinity whereas u ~ g^'^, so a further prescription 
is needed, i.e., 

/ udq""^ (M-n)dg + X(s,A)|,=_i/2, X(s,A)"=2:M {V{q) + Xy^dq, (11) 

Jo Jo Jo 

the latter meaning the analj^ical extension of the formula from its convergence domain 
{Res > 1/A^, A > — inf l^}). As an example, 

+ CXD h 

U{q)dq = ^A^ ioTV(q) = q^. (12) 

4 sm TT/i 

In general, however, X^ develops a pole at s = —1/2 with residue —Z{0)/2, requir- 
ing further regularization in general. All present calculations are carried out under the 
simplifying assumption Z{0) = 0, a special case which still includes whole classes of po- 
tentials unconditionally, for instance: V homogeneous, or purely odd, or purely even with 
N multiple of 4. 

The new normalization (10) is then well defined, and it commutes with spatial trans- 
lations. 



1.4 Basic identities 

Under the above notations and restrictions, very simple identities connect the spectral 
determinants and the naturally normalized solution: 

D-{X)^MO), D+{X)^-^i{0). (13) 

(The proof is an adaptation of the arguments in |Q, Apps. A and D). 

Next, following [||, we continue eq.(P in the complex g-plane down to the rotated 
half-line having the adjacent Stokes direction, namely [0,e~"^/^oo) with 

(fi = — {spectral symmetry angle). (14) 

By simple complex scaling upon q, 

^™=V'e-.A(e'^/^g;tTW) (15) 

provides a recessive solution to eq.(0) in that Stokes direction, where 

expresses an action of the discrete rotation group of order (A^ + 2) on the coefficients; 
equivalently it acts upon the potential V, mapping it to V^'^^, then V^'^\ ■ ■ ■ (now complex 
potentials). We call L the order of the effective symmetry group (L = A^ + 2 generically, 
while L = ^ + 1 for an even V). 

Now, the Wronskian of the two solutions ip,'ip^^'^ of eq.(|l|) is a constant, which can be 
evaluated as g — >• oo because the domains of validity of their respective asymptotic expan- 
sions (H) overlap, and also expressed at g = by means of the respective identities (jT^) 
for the potentials V and V^^^. Matching the two calculations then yields the fundamental 
bilinear identity 

e+'^/^D+(e-'^A,wW)D-(A,w) -e-'^/^D+(A,i;)D-(e-'^A,t;W) =2i. (17) 

We stress that our approach skips any other matching of solutions, such as those in 
connection problems between nonadjacent Stokes directions, which yield nontrivial Stokes 
multipliers... 

We also mention that similar identities have appeared in quantum integrable systems, 
as "quantum Wronskians" [[10 |. 



2 Exact quantization conditions 

For a homogeneous potential, we currently believe that the functional relation ( p!?]) and the 
asymptotic law (H) (imposed to some o(l) accuracy) together suffice to specify the whole 
spectrum exactly: indeed, we can recover the spectrum within each parity sector as the 



fixed point of an (apparently) contractive mapping built upon those data alone [|| • 
Now we start to extend this approach towards more general polynomial potentials, still 
keeping Z(0) = (a broader generahzation is under study). 

Guided by the homogeneous case |^ , and in order to obtain equations for the eigenval- 
ues of a given parity, say even, (the odd case admits a symmetrical separate treatment), 
we set A = —E2k in eq.(|l7D and its partner with v^~'^', then eliminate the D~ , to find 



Z}+(e-i'^A,t7[-i])/^^(e+'^A,t;W)U=-s,, = -e'^/^. (18) 

Whereas the procedure for V{q) = g^closed upon itself at once, this general one invokes 
the complex-rotated potentials V^'^^'^ and, step by step, all L distinct complex potentials 
y[^l as well, their (complex) spectra {E2I} constituting new unknowns (they are just 
analytical continuations of the initial spectrum in complex v-space, but we cannot exploit 
this feature constructively). 

Next, the absolute phases in eq.(p^ (and its cyclic permutations on the -yM) are 
determined by reference to the homogeneous case (see |§|). This operation is essential 
to get Bohr-Sommerfeld-lookalike formulae, governed by an explicit quantum number 
k = 0,1,2, .. .. Our final result, once both parities have been merged again, is then 

logD±(-e-^Ef ,#'"") -log/^^(-e+'^4'l, #+11) = (19) 

for k ^^^^ , i integer mod L 



TTl 



, 1 A^-2 
k + -± 



2 2(A^ + 2) 



the branch of logZ)(A) being defined by continuity from (A = 0, v = 0). 

In each parity sector independently, eqs. ([19D form a system of constraints tying each 
eI at fixed i, k to the two other spectra {£"]• }, {Ej } (whose zeta-regularized prod- 
ucts build the left-hand-side determinants). The system of all such interacting points is 
better displayed, in proper relative positions if each spectrum is first suitably rotated (we 
then call it a 'chain'): the i^^ such chain is then the set {e'^*^ E^^}, i = 0, . . . ,L — 1 (mod 
L) {k being either even or odd, depending on the sector). Each point of the t^ chain 
interacts with every point of the two adjacent {{C. ± 1)**^) chains, through the logarithm 
of their difference entering eq. (|19|) . 

The t^ equation is now a complex one whenever its unknowns E\, are themselves 



complex, i.e., for £ 7^ or L/2: eq . (|T9|) then remains a formally 'complete' system of 
mutual constraints for the unknowns E\^ . As in the homogeneous case we then surmise 



that eqs.(P!9|), which are exact, also constitute genuine quantization conditions, i.e., they 
have the capacity to specify all their unknowns provided the asymptotic condition (0) is 
also enforced upon each spectrum {-E"}.!} separately (using the rotated v^^^). Whenever this 
property can be confirmed, eqs. (^) will then determine the spectral problem analytically. 

We will now report numerical experiments upon quartic potentials which numerically 
confirm this conjecture in some regions of parameter space, by achieving effective com- 
putations of the spectra out of eqs.([T9|). 

(To concretely apply eqs.(|l9D, we evaluate the functions log D^ through the formula (H) 
(adding higher correction terms to accelerate convergence), and we use Newton's method 



to solve eq. ([19D for each El in turn, but only up to some finite fcmax beyond which levels 
are semiclassically assigned once for all; the height of the cutoff fcmax regulates the final 
accuracy.) 

2.1 Even quartic oscillator 

We illustrate the preceding formalism through the quantization of the levels for the even 
potential V{q) = q^ + V2q^ on the real hne. Parity implies: that the spectrum splits 
into even and odd sectors, corresponding to the Neumann and Dirichlet problem on the 
half-line, respectively; and that the order of symmetry is L = 3 (instead of 6 generically) . 

To the initially real V2 giving the real spectrum {E^} are then associated: a com- 
plex spectrum {E^} for the coupling constant jt>2, and {E'l} for j^t>2 (where j = e^'^^/^). 
The three chains {E^}, j{i?^}, ]^{E'l} are shown on fig.l for some values of the har- 
monic coupling constant V2- By reality symmetry, there are only two independent equa- 
tions/unknowns (one real, and one complex). 

The semiclassical information (0) was input up to 6 terms {0{E^^/'^)) to make the 
accuracy grow faster with respect to the truncation fcmax- 

The numerical procedure is to solve eqs.(|l9|) in some cyclic sequence and manner that 
will achieve a contractive-looking iteration, as in the homogeneous case; there are now 
many more implementation details left open, and we tried several. 

The homogeneous case {v2 = 0) enjoyed full ternary rotation symmetry, and also 
yielded strong contraction ratios (< 0.4) |Q, |^. The main change seen when f 2 7^ is that 
certain iteration schemes excite new chain fluctuation modes that are quite asymmetrical 
(even about the fully symmetrical f 2 = chain) and much less stable (contraction ratios 
~ ±0.9). This holds for the naive (and most symmetrical) deformation of the f 2 = 
iteration, which at each step addresses each chain in turn but refreshes all of their values 
only after a complete £-cycle. Those 'bad' modes ruin the quality of convergence already 
very close to 112 = 0. 

By contrast, we empirically found that an iteration scheme alternating between the 
real and one complex chain and immediately refreshing their values at each step, remains 
fairly stable near f 2 = while appearing to converge towards the correct spectra. It 
nevertheless deteriorates as V2 grows, becoming totally erratic above V2 ~ +2; hence it 
cannot be carried towards the harmonic limit. By contrast, the scheme looks much more 
robust in the negative direction (double- well region), showing no perceptible degradation 
down to V2 ~ —10 ! Its contraction ratios are plotted on fig. 2, left (for ^2 < 2 they were 
estimated from the iteration results and by the diagonalization of the linearized dynamics 
near the fixed point, and for f 2 > 2 only by the latter method). 

Some explanation of the results comes from fig.l: when V2 grows, the complex chain 
and its conjugate become almost degenerate towards low quantum numbers (a sort of 
complex tunneling effect), increasing the logarithmic pair interactions and the instabili- 
ties (the linearized-dynamics matrix entries grow even faster, like [mutual distance]"^). 
An obvious improvement is then to decouple the interaction between any two complex- 
conjugate chain points by enforcing this symmetry as an explicit constraint (immediate 

6 



refreshing also has to be kept); this achieves a much more uniform hnear stabihty indeed 
(fig.2, right). However, global instabilities still arise, albeit more slowly, and we still lost 
all convergence beyond V2 ~ +5 (those instabilities primarily hit the Newton root-finding 
algorithm, so there is room left for remedies). 

The preceding numbers referred to iterations in the even-parity sector. As fig.2 shows, 
the odd sector tends to behave more stably. A conceivably more robust scheme might 
then be to confine iterations to the odd spectrum, then to solve for the even determinants 
from the coupling relations (|l^ instead. 



2.2 Extension to non-even potentials 

In we showed that the exact quantization formalism was fully valid for homogeneous 
potentials of any odd degree iV, and specially for A^ = 1 which has Airy functions both 
as solutions and as spectral determinants. The paradox in the very regular behavior of 
this case is that the underlying potential V{q) = \q\ is not even once differentiable at 
g = 0. This exemplifies that the 'even' and 'odd' spectrum distinguished by the exact 
quantization mechanism are just the Neumann and Dirichlet spectra. Parity properties 
of the polynomial defining V are then irrelevant: the underlying even potential can and 
will be ^(Iq'I) in all cases. 

However, a generic non-even polynomial V{q) has Z{0) 7^ 0, a case which our approach 



still handles incorrectly. To test the exact quantization conditions (|T^) upon quartic non- 
even polynomials, we therefore tried a few with Z{0) = 0, using 

Z(0)[=(2.)-5„]=-| + ^-| for AT = 4; (20) 

for instance, V{q) = g^ — g^ -|- |g|/8 yielded results comparable to the previous ones for 
even systems. A further validation is also brought by the next computation. 

3 Exact wave-function resolution 

That extension of the exact quantization formalism beyond even potentials is important 
for the following application. If no parity property is imposed upon the polynomials V{q) 
with respect to the endpoint g = 0, then this endpoint can be taken to an arbitrary value a 
(on the real line, at least), thus restoring translational invariance effectively. Equivalently, 
is kept as endpoint but the potential gets shifted to V^(g) = V{q + a) — V{a) on the 
half-line (becoming V{\q\ +a) — V{a) on the whole line). It further happens, in the special 
case ^(g) = g^ + f2g^, that our present restriction Z{0) = remains fulfilled by all the 
potentials V^(g) (simply using eq.(^)). Hence this provides another opportunity to try 
the current formalism. 

So, we now use eq.(ll4D at g = a instead of and from right to left: then it states that 



the naturally normalized solution of eq.(|I|) at g = a, resp. its first derivative, equal 

ip{a) = D'{V{a)+X), ip'{a) = -D^{V{a)+X) under the potential K(g) on [0, +00). 

(21) 

In turn, these determinants are canonically specified as the zeta-regularized products 

over their chains of zeros. And these are the eigenvalues of V^^^{q), hence they obey 



exact quantization conditions of the form ([T9|) and, hopefully, this makes them reachable 



by some convergent iteration scheme. In this sense, eq.(pTD will provide an algorithmic 
specification to solve eq.(|l]) (for arbitrary A). 

We illustrate this program with a calculation of the ground state eigenfunction for the 
homogeneous case V{q) = q^. The shifted potential is then K(g) = q^+Aaq^-\-Qa?q^-\-Aa^q, 
and 'ip{a) is the value of its determinant D~ at the point a^ — Eq, where the eigenvalue 
£^0 ~ 1.06036209 is part of the input here (but the levels themselves are the output of a 
similar exact calculation using V{q) = q"^). 

The exact quantization conditions for Va{q) now involve L = 6 distinct chains (except 
at a = 0, where L = 3 by parity symmetry): i = 0, 3 are real, (1, 5), (2, 4) being doublets 
of complex conjugates. We only tested ipi^a), for a few values of a > 0. We found this 
iteration sequence: £ = {0, 2, 3, 1} (then cyclically continued, with immediate refreshment 
of each chain and straight enforcement of 1-5 and 2-4 symmetry) to converge indeed for 
a < 1.7 (contraction ratio per cycle < 0.67, better than other sequences for no clear 
reason). Just as above, instabilities suppressed convergence for higher a. Our output 
points are plotted on fig. 3, against the curve produced by a standard integration routine 
and upon which only the global normalization was fitted: the results show an overall 
4-5- digit agreement. 

3.1 Concluding remarks 

We have reduced the resolution of polynomial Schrodinger equations (|I|) (for their eigen/values/vectors) 
to that of a discrete system of selfconsistent exact quantization conditions, eq. (pJ]) , having 
as unknowns (A^+2) countable chains of points subject to asymptotic boundary conditions 
supplied by conventional Bohr-Sommerfeld formulae (^ (still provided the zeta function 
Z{s) of eq.(|D vanishes at the origin.) 

We are also getting numerical evidence that the so reduced problem is effectively 
solvable in some regions of parameter space through iterations which seem to converge 
geometrically. This suggests that the solution is being described as a fixed point for some 
explicit action of eqs. ([I9|) upon infinite but asymptotically constrained chains, appearing 
to be contractive and also to admit robust finite-dimensional approximations. 

On the other hand, since eqs.(|l9D are very tied to the values of the degree A^ and 
symmetry order L, they may be ill-suited to transitional regions where these numbers 
jump (e.g, f 2 — ^ or + 00 in V{q) = q^ + V2q^)- Moreover, our available evidence is 
currently confined to even quartic potentials close to q^ and (from earlier work [Q) to 
purely homogeneous potentials for their spectrum only (but up to quite high degrees). 
Extensions to higher A^, 2'(0) 7^ 0, arbitrary complex q and v etc., and to more general 
differential systems, are all conceivable but they remain open issues. 



Finally we argue that, while the integration of eq.(|l|) by quadratures alone may still 
be a remote goal, another relevant question is how much the set of admissible integration 
methods needs to be enlarged to reach that purpose. We are then finding that zeta- 
regularized infinite products (of order < 1) plus the solving of one type of fixed-point 
equations (which seem to have contractive and other nice properties) appear to be well 
adapted, and sufficient to the task already in a few cases. 

Acknowledgment: we are grateful to R. Guida (Saclay) for putting at our disposal 
his computer program to calculate hundreds of anharmonic oscillator levels with high 
accuracy. 
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Figure captions 

Fig. 1. Interacting cliains in tlieir equilibrium positions yielding the exact odd spectra 
of y'^1 (points marked with the corresponding £-value), for potentials V{q) = g^ + ^2^^- 
As V2 -^ — oo, the complex chains tend to shadow the resonance spectrum of the potential 
— \v2\q'^', as V2 -^ +00, the complex chains tend to shadow what would be the i = 1 chain 
for the potential V2q'^- 

Fig. 2. Contraction ratios of iteration schemes converging to the spectra of V{q) = 
q^ + V2q'^ (+: even spectrum; o: odd spectrum). Left: numerical estimates for simple 
iteration scheme with immediate refreshing; we note that the numerical values never reach 
exactly unity, even though no definite conclusion can be inferred (no error estimate at all 
is implied in our data, but 2-3 digit stability is typical). Right: numerically estimated 
moduli of same ratios but for an iteration scheme further decoupling conjugate pairs (here, 
different parts of the curves may correspond to different branches of linearized eigenvalues, 
some being negative or in complex pairs). 

Fig. 3. Naturally normalized ground-state eigenfunction for the homogeneous quartic 
potential g^. Curve: computer integration of the Schrodinger equation using the NAG 
routine D02KEF; its rescaling factor was the only number fixed by a fit. +: points 
obtained by iterative exact quantization; o: numerical estimates for the corresponding 
contraction ratios. 
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